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preface 

This mejmo ts aat Intended to De of tiitorlal value. Xt 
presumes t^ general kn^o^^Xed^e or POOTRA-N., suon as can be 
obtained t^cjus A FOkHiMH IV PRIMER by E. I. OrganlcJc. 

Definitions and notation 

For clarttj, e^tamples of Faerj?2IAN eoxxrce st^iteiner^ts Mill 
be written m uppey c&ae, 

pyoce<3i\kre losano fanctlen or mibroutlne^ (External proca- 
durea are aoraatimes called subprosrame . ) A function retui?iis 
a value, and is called by a reference in an arithsiefctc 
erp3?93»ion. A subrcfutlne aoe» not reUirn a value # and Is 
called by ttee GALL statener^t. 

Mode {also knomi as type) is one of Integer^ reari, 
logical, cojsplex^ &r}id. dajblepreclsion. With very few excep- 
tional all data processed by FORTHAN programs are tn otise of 
these modes c (One exception is a procedure naioe arg\sment to 
a proced\2re,>) The iRode of a constant is determined by Ita 
formo The mode of a variable Is determined by explicit mode 
declarations^ Implicit inode declarations^ or the IJKLMN 
rule^ In that order. -fhe mode of the value returned by a 
function Is determined In the same vmy^ except for intrinsic 
and library fwjctlons. 

The XMLm rule is that variables and functions whose 

names begin with X^ J,, K# h, M^ or K are integer mode^ and 

others are real,. This rule jnay be changed by an IMPLICIT 
statement* 

A dimension declaration is a declaration consisting of 
the array name followed by a llst^ in parentheses^ of the 
iBaxlmum values of each s^jibscrlpti It may appear in 
DIP^NSION,, LOaiGAL^ REALi, DOUBLBPHECISXOH, INTEGER* (X)MPLEX^ 
and comoii statements o 

A data declaration Is a declaration specifying the 
initial contents of a variable or array^ and may appear in 
DATA, DIMENSION, LOGICAL, EEAL^ DOUBLEPHECISION, lOTEGERij 
and COMFLSX statements » 

An ASF is an arithmetic statement function (internal 
function defined in a single statement which looks like an 
arithmetic assignment statement with subscripted left side)^ 



Sow^ee program fornsa" 

fi »t;aten*ent la a string of character's 'ap to a carrlEge 
retvirn, not counting carriage return a hidden by center diot or 
in character strings. A centerdot causea all characters up 
to and including the next cax^riage return to be ignored^ and 
Is us«d for continuing statements on more than one line* 
There le no limit to the mimlaer of ccntinuatlon lines tlmt a 
statement way in^ve. An empty statssjent or statement with Cj 
period, or slash as the first character is ignored. Center- 
dot may not ice used to contim;ie consaents. Except in 
character strings, space and tab are ignored. There are no 
other restrictions on fornisit or line position » Specif icallya 
statement numbers and the rest of the statement do not have 
to occupy specific ^fields^. Statements nay begin at the 
left inargin with no tab, but beware of CONTINUE, COMMON^ 
COMPLEX, Gft-LL, and other stateaa^fnts beginning with C- These 
must be preceded by a space cr tab if they have no statement 
number, to prevent their being treated as coinments. 

Letters may be in upper or lower case* Brackets nsay be 
used in place of parentheses* The case of letters is 
ignored, i*e«, the names ABC, ABc, and abc are one and th<© 
3a.me« 

The first non-enipty line of a program is the title. It Is 
printed on the typewriter at the beginning of compilation 
and is not a statement* 

Character strings aay be written with '^numberK^ notation 
or asay be enclosed in single quotes. In the former case,.* 
case shifts are counted as characters ^ and they are not 
filtered. In the latter case, case shifts are filtered ,« and 
the string is assumed to begin and end in lower casec Single 
quotes nay be placed in the second kind of string by using a 
pair of consecutive single quotes. All characters are legal 
m either kind of string. Strings in other than format 
statements are left Justified if they contain more than 3 
characters • They are right Justified and considered to be 
Integer constants If they contain 3 or fewer characters. 



Arithmetic expareealona 

A symboX Is a string of xjp to 6 letters or aiglts^ the 
first of which 18 a letter o A syaboi oay be used as a 
varial>l^, array n&ciej function naiae, subroutiwe narR©^ CCW^ON 
block name J or KAMELIST name. Except for COM«ON block names, 
a symbol nay be only one of the above ^^nlesa It Is a 4uBsny 
cymbal for an internal psrocediire^ in which case it may have 
one laeanlng within the procedure and another meaning else- 
where « 

A variable Is a nan^ for a reglatei* in which on^ data 
itea i»y be stored during execution o The mode of the 
variable determines the type of aata that laay be etored. 

*n luteger dat^an 1« a signed integer trith laaagnltv^e not 
greater than I31071 decimal or 377777 octal. 

A real dat^si l« an approxljoatlon to a real nuaiber, it has 
nagnitude between lOt--75 and lOtTO^ or is zero^ 26 
bits (about 7 1/t decimal digits) of precision are 
Icepto It is stored in two consecuatlve words o 

A doubleprecislon datum has isagnlttide between iOt«39000 

and I0t39000^ or is sero. 69 bits (about 2% decimal 

digits) of accuracy are kept. It is stored in 5 
consecutive wox^ds,, 

A logical datum has the value .TRUE. (stored as 1} or 
« FALSE ^ (stored as 0)0 

A complex datvim is a pair of real data^ representing the 
real and imaginary parts of a complex number. It Is 
stoi^ed In four consecutive words ^ 

The mode of a variable is determined by an explicit mc?de 
declaration (IHTEOER,, REAL, LOGlCftLi, COMPLEX, or 
DO0BLEPRECISION statements}^ or, if none is present ^ from 
its first letter, using an IMPLICIT declaration or the 
IJKLMN rule. 

An array is a set of variables named by a single symbol c 
Individual items are selected by a string of subscripts. 
Each array must appear in a dimension declaration^ giving 
the ntaaber of subscripts and the maximum value of each* Any 
reasonable number of subscripts may be usedo Each subscript 
my vary between i and the limit given in the dimension 
declarafcion. The items of an array (each of which is i^ 2. 
^p or 5 words long depending on the mode) arc stored in 
ascending locations in memory* the first subscript "varying:" 
most rapidly.. 

example « DIMENSION v7J(g^3) 

The elements of the array^ in order of increasing 
address j> are 

jjCiAl jjC2,i) jjCi^g) ^jCt,^} 



An array reference is an array name foXlowcd by a 
parenthesized list of subscripts. Subscripts may be of 
arbitrary complexity * The number of subscripts irsust match 
the ntmiber of dimensions given in the declaration^ or roust 
be one. In the latter case, it is a "linear subscript** and 
selects an item from the GiVT&y as a linear list* 

The mode of an array is the mode of each item in the 
array and is determined in the same ii^ay as the mode of a 
variable « 

Constants 

An integer constant is an optionally signed number^ octal 
number, or character string of 3 or fewer characters* An 
octal number Is a nuaiber followed by the letter 0, A 
character string is enclosed in single quotes or is preceded 
by a number^, giving the length of the string, and the letter 
H. The characters are right Justified and filled with zer<s 
on the left. 

exai^ples - 123 -iR3o *iha '^foo*' 

The letter in an octal constant may be in front of the 
digits in a data declaration, where the constant cannot be 
confused with a symbol. 

A real constant is a number with a decimal pointy an 
exponent preceded by E, or both^ The exponent may be signed, 
as may the entire constants If there Is no e:s:ponent* the 
number of digits must be ^ 9 to distinguish it from a 
doublepreclslon constants 

A doublepreclslon constant is a niaasber with a decimal 
point, an exponent preceded by the letter Dj> or both* The 
exponent may be signed, as may the entire constants if there 
IB no exponent, the number of digits must be > 9* 

A logical constant Is oTRUB. or iB for truth, .^LSEr or 
OB for falsehood 6 the forms T and F are also permitted in a 
data declaration, where the constant cannot be confused with 
a symbol « 

A complex constant is a pair of real constants separated 
by a comma and enclosed in parentheses « Each part may be 
signed, as may the entire constant* 

Note — The form left par en, real part, comma, imaginary 
part, right paren, is legal only for constants « |a,B+C) is 
illegal. Use CMPLX(A,B4-C) Instead « 

A long character string constant is a string of vaore than 
3 characters. It Is not one of the 5 data types, and is 
legal only in data declarations, procedure arguments, and lO 
lists In a dat^ declaration, it is treated as a list of 



IntegeTs. in a pr'oeed^Jre ayg\a5Q^r»t it ia ts^eatisd as an ar.?*si>% 
A long character starring Is lef^ Justified and filled wlUh 
iseroes on the right. 

A etatemsnt ia^el constant is a stateraent numbcsr foiXoirea 
by the let^ter S. It is not on$ of feh« 5 datn. tjpeta^ and l8 
lagal only In data declarations and subroutine apgumsrsts. In 
a data declaration It 1» used to inltia-llse a variable to be 
iJiW2d In an aaalgnefi GO TO« In * sul&routlsc* aipgtainent it is 
ueed to specify an a^are^a for an errcvr retmam^ The letter S 
may precede %h& number in a data declaration j, inhere it 
cannot be con/used with a aj/nbol^ 

A functlc^ reference ia a functlcai name followed by a 
p3u:*entheBlzed argumen*; liat., The laode of a function refe3?« 
ence is the awde of the function naBWi which^, accept for 
intrinsic fiarictiona and standard library function** 1^ 
detemlned in the saai3 way as the m.>dte of a varlabliie 

An expression is a constanti variablef, array reference^ 
function referencei> expression enclosed in parentheses^ 
expression preceded by a miary operator^ or expressions 
8ei;^rated by binary operators* 



IJbe arlthwetlc operatsors apre listed below 
operator(3} acceptable sic3de(Bj 






I* 








logical and 


L 








logical or 


1, 








in^iee Cb or not A) 


L 








logical exclusive get 


I 








logical eq'ijiivalQiiee 


h 








logical negation 


I 








bitwise and 


I 








bitwise GP 


I 








bitwise ©KcltiBlve or 


I 








ones coinpX^itent 


X 








renalrxl^r 


2 








X shifted left by -f 


2 








X shifted right toy ^ 










X 3:*ot&t«d left by S' 


I 








1 rotated right by •! 


I:. 


a. 


D* 


C 




X. 


a. 


D, 


G 




I^ 


E* 


B* 


C 




X, 


H^ 


D, 


C 




I* 


R. 


D^ 


c 


exponentiate 


I. 


R* 


D 




absolute value 


X. 


B. 


D 






I. 


E^ 


D 






I, 


H. 


D* 


c 




X. 


R* 


B, 


c 




I. 


a^ 


D 






I. 


B, 


D 







- C^^najpy or binary) 

X 

/ 

^ABS, (xinary) 

«<Xt]&4 ^ 

Where parentheses are not %xa^ to direct the order of 
computation, operations are performed in the order shown 
below „ 






ABSo unary - 

XaSe. oX^o «HSd uHH. 

K. 

A 

oEV, "• 
.Vo V 

/ .REM, 

X 

+ binary - 

eOTo oO, > oGEc ^ .EQ. 

<,AND« 

,0R, 

•THENo 



lone first) 



.E. 



JN£ 



XT 



O C<*i< 



(done lest" 



Where two or mere non-comaautatlve operators oV the eame 
prlorlt;^ levtl are U3ed« the operations are perfcyrmeai jfrom 
le^t to Tt^t except for exponantiatioiii i«hich is performecl 
from right to left. 

All arguments to an operator are first converted to an 
acoeptftble 8K?de t?y the snallest possible subset of the 
follovlng eonyerslonao 

doublepreclslon ^ re$.l 

H^^ leea alg^^lflaant bits of the fraction eu^e dlscaiH^ad, 

real -»- (ioublepreclslon 

Zeros are appended to the fraction « 

doubleprecialon ^ integer 

1!t)e number la truncated toimrd aero^ It is not rounded, 

integer •*• doublepreclsion 
The fraction is zero. 

coroplex -* real 

The real part ia used^ 

real •♦ complex 

Zero is used for the imaginary part^ 

real -* integer 

The number is trtmcated toiiiard zero^, It is not rounded c. 
2.999 -* ^ '«2.999 ^ -2 

integer -> real 

The fraction is zero^ 

integer •* logical 

Zero is converted to oFALSE* Anything else is converted 
to ^TRUE« 

logical •♦ integer 

oBSJlLSE, becomes 0^ .TRUE, becomes io 

Except for exponentiation with Integer second argument « 
all arguments are then converted to a consistent mode 
according to the following 5 steps. 

If there are any doublepreclsion arguments/, all integer 
arguments are converted to doublepreclsion „ 

If there are any real or coo^jlex arguments, all integer 
arguments are converted to real^ 

If there are any complex arguments* all doublepreclsion 
arguments are converted to real,, 

If there are any doubleprecision argmaents^, all real 



argy©fent5 t-re eonvorted to doubl©iareci»lcni. 

If there are aarj compl.&x. argiirttenfes, aDll real argument* 
are conv^^pted to coi»piQx> 

Whenevei* ^wo or more argiuroente tutderso th« saiBd conver-* 
alon, the arithmetic opera tlcm le performed on them first, 
and the seeault convex^ted » For ex^ungple ah^I-^&^.^-km-K is 
rearranged as A4B+OFL0ft.T< X-hI-^^IC) e The intege;? argianenta are 
added a a integers and the^i converted « 

When the seccmd argument for exponentiation la integer i, 
It is not converted to ne&lce It consletent ulth the flrat^ 
that 1&|» the operations doublepreclolont Integer^ 
realrintegerj and complextlnteger are performed without 
conrereion, 

Except for connectives (•OT^j !,EQ.# etc.) and ejqsonen- 
tlatlon, the mode of the result la the mode of the 
arguments « The result of a connective Is alimys logical » The 
result of exponentiation has the mode of the first arguments 

Except In 10 lists and arguments to procedures » expres<« 
slons are alimys converted to the required mode« For 
example^ the right side of an arlthn^tlc assignment state- 
ment Is converted to the mode of the left slde^^ a computed 
00 TO Index Is converted to Integer, and the test expression 
In a logical IF is converted to logical ^ There Is no *mlxed 
mode" error* 

When the result of an Integer operation Is not Integral 
(this may happen In division or exponentiation with negative 
exponent) it is truncated toward zero, in the forms 
realtx*eal and douhleprecislontdoublepreclslonp the base must 
be non-negative « OtO Is undefined in all modes c 



Executable statdments 

Statement labels (atateisent nunibers} 

Any executable statesoent may be preceded by a label > 
which la a decimal integer between 1 and 99999 > Statement 
labels are used to generate addresaee for GO TO and IF 
statements^ and to Indicate ends of DO loops « labels are 
ignored on non-*executable statements other than FORMAT, 
Labels on F0IU4rv7 statements must not conflict with labels on 
executable statements. 

Arithmetic assignment 

The arlthiaetlc assignxaent statement consists of a vari- 
able name or arztiy reference^ an equals sign^ and an 
expression* The expression is evaluated and stored in the 
indicated variable or array element. If the left side is an 
array reference, its subscripts may be of arbitrary com- 
plexity and follow the rules for araray references in 
expressions. The expression to the right of the equals sign 
is converted to the mode of the variable or array oa the 
left by the smallest possible subset of the conversions 
given in the section on expressions » 

Simple 00 TO 

The statement GO TO followed by a number causes transfer 
of control to the numbered statement • 

Computed GO TO 

The statement GO TO followed by a parenthesized list of 
statement numbers (separated by commas) followed by a comma 
followed by an expression* is a coinputed 00 TO, Control 
transfers to the statement number whose position in the list 
is the value of the expression. The expression may be of 
arbitrary complexity. If It is not integer ^^ it is converted 
to integer. If the value of the expression is 1, control is 
transferred to the first statement in the list, etc. 

example « GO TO (10,30»25)* J+K/2 



ASSIGN 

The statement ASSION^ etatement numbert T0« variable cr 
array reference » storea a pointer to the specified statetoent 
in the variable or array element. The variable or array 
shoiild be integer* This statement is used In conjunction 
mth the assigned 00 TO statensento 

exaaqpae -• ASSIGN 4 TO jr(K+L} 

Assigned ao TO 

The statement 00 TO followed by a symbol or array 
reference f ollcired by a c<»nna followed by a parenthesized 
list of statement numbers > causes transfer of control to the 
statement whose pointer was last stored in the variable or 
array element by an ASSIGN statements The parenthesized list 
of statement numbers is required and must contain all 
possible statements to which control may be transferred « The 
conipiler uses thes information internally » 

example - 00 TO J{N+i), (3*^*5*14) 

The use of an array for assigned 00 TO^s facilitates certain 
"pushdown stack" operations. 

Arithmetic IP 

The statement IP followed by a parenthesized expression 
followed by 3 statement numbers separated from each other 
(but not from the expression) by comn^s, is an arithmetic 
IP, The expression is converted to integers real, or 
doubleprecislon, and tested. Control transfers to the firsts 
second or third statement if the value is strictly negative;^ 
zero, or strictly positive, respectively. 

exan^le « IF (j-a) %0,1L3,10 

Logical IP 

The statement IP followed by a parenthesized expression 
followed by kny executable statement other than DO^ is a 
logical iPo The expression is converted to logical « The 
statement following the expression is executed if the value 
is .TRUE. 

exanqple »• IP (P ,GR. Q) J « J+1 



DO 

A DO loop Is h^gvm by the statement D0» a statement 
namber* symbol, equals sign, and a or 3 expa?easions sepa- 
rated by QQmm&. The loop ends Just after the tormSml 
©tateaent^ whose number is given in ths DO statement, Th^ 
symbol before the equals aign in the BO statement Is the 
indexe The exparessloijs after the eqmls sign are the initial 
vaJue, test value^ and incrementi, respectively^ They may be 
of arbitnary complexity. if the Increiaent is misalng. i 
t integer) is used. Vhan the DO statement is execiuted^ the 
index la set to the initial value j» and subsequent statements 
up to axid including the fcermlnal at^tement (the '^range*) are 
execmted repeatedly. The Increment is added to the ina0K 
after each iteration. If the result is less tl^n or eq^l to 
the test values the range is executed again ^ Otherwise it 
proceeds to the statement after the terminal statement^ Th# 
range is always executed at least once. Mode conversions are 
performed where necessary to do the addition and testingo 
The testing is done as if by a .is, connective ^ so the 
argents are co?iverted to real^ integer f, or 
doubleprecisiono Because of the i^y the test is performed. 
fcht^ increment should alleys be positive. The terminal 
sta tenant may be any executv^ble statement other than 00 TO? 
STOP^ RETlfRH fisnless wider logical w}s or another DO. It is 
permissible to Ji^mp in and out of DO loops at any time. The 
last value of th€ index (after being incremented) is 
preservea upon nornial exit from a do loop. DO loops my be 
nested to any reasonable depth « 



Ithtf COmWOK stateisient ao^s nothing. It is u&ed aar the 
textdnftX stateiMnt of a DO loop In order to be able to go to 
thB end of the range of the Xoop<, 

BftTISE ajRd STOP 

The atateoiont i^VBK or STO?» folXowea ^ nothing^ an 
octal niintoer, cxr a character string enclosed in aln^Ie quote 
mtrkSj eaueee execution to atop» If the statement has a 
ntii*>er ov etring. It Is tjrped out flrat. The nwber, if 
present, will be printed In octal. 7he STOP statesient 
termlnatea execution with the dam Instruction^ The PAUSE 
etatenent tdalts for any character to be typed In and then 
proceeds » If a character string is given, it must be 
enelQsed In quote markSo •tf^ foroat la Bofc permitted {tt 
would be ambiguous). The aB>w» effect as a STOP statement aaay 
be achieved by running Into an END statement (running off 
the end of the program or procedure) or by an end of file on 
a READ statement when no end-of-file address Is given. 

examples •- PAUSE PAUSE 777 PAUSE ^foo^ 

CALL 

The statement CALL follwed by an Internal or external 
subroutine name, followed by an optional parenthesized 
argument list, calls the subroutine. The form of the 
argua«nts is discussed in detail in the section on proce- 
dures, 

RETURN 

The statement RETURN followed by nothings, a number (for 
non«standard subroutine retvima)^ or a parenthesized expres- 
sion (for function returns) «> returns ft»om the innermost 
procedure. This statement is discussed in detail in the 
section on procedures. 

READ and WHITE 

These statements are used to effect transmission of data 
to or from the outside worlds READ or WRITB is followed by a 
parenthesized list of one or two items^ followed by the 10 
list. No coffioa precedes the 10 list^ The first item in 
parentheses is the \init number^ which is an expression of 
arbitrary ccanplexity. It is converted to integer if neces- 
sary. The second item, if present ^ specifies a PORWAT 
statement number or ari«y to be used for ^variable format"^ 
These statements are treated in detail in the section on 
input/output « 



Declarations 

Th« modi© declarations (RBAL^ LOaxCAl., COiySFLEX, and 
BOUBLEPRECISION) 

Thece declax^tlone are used to explicitly apeciry fche 
mode of a varl&t>le» duxmy syratool, or function name» The mode 
word la followed lay the namee to be declared « which are 
separated from eaoh other j> but not from the mode ward, by 
comroaB. Mode declarations uay be placed anywhere^ except for 
declarations on duunmy symbols for internal procedures, which 
must be Inside the procedure definition. 

If a function is called in a program, a mode dec2.aration 
must be given unless it is a library or intrinsic function 
or the correct mode is given by an IMPLICIT declaration or 
the XJIOMM rule. In a FUNCSTION program which defines an 
external function ;. the mode must be declared « either with an 
explicit mode declaration or In the PUWCTION statement 
itself, unless the correct mode ta given by an IMPLICIT 
declaration or the XJTKLMN rule, in a program in which an 
internal function is defined and used^ its mode must be 
declared either inside or outside of the definitions^ or in 
the INTERML FJNCTION statement itself^ unless the correct 
mode is given by an IMPLICIT declaration or the IJKLMN rule. 

A mode declaration may give dimension Infornation for 
arrays by replacirig the symbol name with a dimension 
declaration (following the symbol name by a parenthesized 
list giving the maxlmaam value of each subSGrlpt) . 

example - LOaiCAL A, bC4»5,£}, C 

A mode declaration rmy specify initial contents of 
variables or arrays by following the symbol name {or 
dimension declaration) by a mt& declaration, which la 
enclosed in slashes^ 

example « LOQK^AL A/TAB(4,f^^g)/40x.PALSE.Ac 

The DIMENSION declaration 

The DIMENSION declaration specifies the size arid sub-^ 
script function for an array, it consists of the word 
DIMENSION followed by a list of dimension declarations, 
separated from each other by comnaas. Each symbol to be 
declared is followed by a parenthesized list giving the 
jcsaximum values of each subscript. Subscripts way vary 
between 1 (not o) and the given limit. Inclusively. Except 
in dimension declarations for dummy symbol s^, the limits must 
be integer constants « since the compiler needs to know how 
much space to resei^e. In dimension declarations for dummy 
symbols* the limits may be variabilis <. This feature^ called 
"adjustable dimensions** ^ facilitates coding procedures that 
manipulate arrays of different siz^s. The limits in the 
dimension declaration are dumny symbol©, and the corre- 



spending argumenfcB aro nuonbers by which the calling program 
specifies the slz© of the array toeing transmitted* 

ex&mpl0 ~ subroutine to transpose a square tnatrix 

SUBROOTINE I'R/INCA, SI2E) 

INTBGEH SIZE 

DIMENSION MSXZE, SIZE] 

DO 13 I « S, SIZE 

DO 13 J « 1* I-^ 

K «A{I, J)^ 

A{I, J) • kUs I) 
2.3 A{J, 1} « K 

HETURN 

When this subroutine Is called with a square matrix and 
length of a side as arguments, the subscript function and 
terminal value for the outer DO loop will be adjusted to the 
size of the array,, 

warning - the subscript limits (except possibly the last) 
in the dimension declaration for a dtjanray array neist always 
match the corresponding limits in the original array, 
calling the above subroutine with a second argument smaller 
than the actual size of the array In an atten^t to transpose 
Just the upper left corner will not workc 

A DIMENSION statement may specify initial contents of an 
array by following the dimension declaration by a data 
declaration* which is enclosed in slashes. 

example - DIMENSION ?00(5}/As .2^ .3# A, 3/ 



The COMMON declaration 

A CQIf^N declaration specifies that c^rtalrs variables and 
aritiya are to be placed in named tolocks of coramon stcsrage, 
where they can be linked by the lc«3id$r among several 
programs » 

The wocrd COMMON Is followed by symbols and block names^ 
the latter enclosed In slashes « The block names are looked 
at only by the linking loader;, and imy conflict with other 
symbols « They must conform to the requirements for PORTBftN 
symbol So 

The symbols to be placed in common storage are separated 
from each other, but not from the word COmON or the 
preceding block name, by commaa, A coraroa separatJlng a block 
name from the preceding symbol Is optional, 

example -- COMMON /POO/ A,,B, /feAR/ C.aD /mm/ E,l? 

Each symbol Is a»slgn«^d to the block whose nem© most 
recently precedes it^ or^ If no block name precedes It in, 
the statement for the preceding block nam© Is null)^ It Is 
assl^rjed to a block called *blank common*^. Within each 
block, variables and arrays are assigned to increasing 
addresses in the order in which they appear <, 

COMMON statements may be placed anywhere. Dummy symbols 
smy not be assigned to common storage, A C0^2M0N statement 
may give dimerision information for arrays hy replacing the 
symbol name by a diniension declaration c 



The DAm declaration 

The DATA declaration apeclfies initial contents of vari- 
ables and arrftye* It aay be placed any^ihere^ The word DATA 
18 folloimd by pairs of variables lists and constants llstSj, 
the pairs being separated from each other hy comn^So In each 
palr^ the variables list is a list of variables^ arrays^ 
array elements, or iiaplled DO**b> separated from each other 
by commas e The variables list nay not contain duasay symbols^ 
The constants list is enclosed in slashes and consists of 
cc«5Stants, or constants with repetition counts^ separated by 
ooiaofts«> Constants must msttch the mode of the variables to 
which they are assigned* 

exainples - DATft A/3.0/,B.C,dA.0^5.0,6.0/,S,F,Q/3>^.0/ 
DATA (Q{l),I-i/eO}AO)d,0,lOx-icO/ 

An intplled DO consists of a l^ft parenthesis^ a list of 
variables, arrays^ array elements ^ or Implied DO<^S|, a cozacna, 
an index rvame (which is always considered integer and is 
unrelated to any other symbol},, eqimls sign^ initial value* 
coflwia^ final value 9 cc^mna, increment » and right parenthesis ^ 
The increment and preceding ccaana may be emitted, in which 
case the increment is one. The compiler Iterates over the 
lists letting the Index vary in the same way as it would in 
a DO statements Implied DO^s may be nested to any reasonable 
depth* Implied DO parameters ami array elwuent subscripts 
are integer expressions that may be of arbitrary complexity* 
They may not contain references to intrinsic or other 
functions, and the only variables permitted are indices of 
implied DO^s contx'olling the expression o 

The constants liat consists of constants separated by 
coamaSo The constants may be preceded by a positlva integer 
and a multiplication sign>j to Indlisate that the constant is 
to be repeated. Several forms are legal for constants In 
DATA statements but not in expressions because they would be 
ambiguous or generate ijiore than one data Item^ 

Long character string constants are legale They are 
treated as lists of Integer constants, one for each 
three characters* 

T and F may be used for logical constants « 

Statement numbers^, in either of the forms i^3S or Sl23^ 
nay be used fog* initializing variables for assigned 00 
TO^s, 

Octal constants may be ^^^Itten as o?7* 4077^ -oTi^ o+*H^ 
or 0-77 » 

For each variables list/constants list palrf, the compiler 
scans the two lists in parallel |> naklng the indicated 
assignments o When either list runs out? the reaalnder of the 
other list is ignored o 



The EXmnmh declaration 

Thfii EXTERNAL etafcement declares symholB that appear as 
arsumerit;8 to procedures as entry points to e^cternal proce-' 
dure* J, to prevent the compiler from treating fchem ae 
variables. This dieclaraiicn is not needed for symbols that 
are already recognizable as entry points. The word EXTERNAX* 
is followed by the symbols to be declared, separated from 
each other by commas, Tho EXTERWAL statement may be placed 
anywhere o 

The EQUIVALENCE declaration 

The EQ0IVALENCS statement specifies that certain vari« 
ables and array elements are to occupy the same locations in 
smmstry^ It may be vsed to place symbols in specified areas 
of coRttDon fitoragej, to give special names to certain elements 
of an array, to patch programs with misspelled symbol s» to 
give a variable different nanvas In different modes* etc. The 
word EQUIVALENCE is followed by a list of declarations* each 
enclosed in parentheses and separated from each •sjther by 
comrsas. Each declaration conr^lns two or tmr^ symbols^ 
arrays, or array elements tmt are to toe assigned to the 
same location. Subscripts for array elements mist be 
constants. An array nan^e is tr«?ated a© fche first element of 
the array, 

example ^ EQUIVALENCE CA,B^C),(DC4}^EC3,t))^CF(E.4)^a) 

This declaration nsay be placed anywhere. It may not be 
used on dufisny symbol©. 

Th?? XKiPLIcrr declaration 

The ISKLm rule may be changed by the XHPJLIGIT statement. 
The word IMPLICIT Is followed by one or more declarations,, 
separated from each other by commas. Each declaration ia a 
mode word tiNTEOER. REftI*, I.OQIQAL,, COMPLEX, or 
DOTBLEFRECISION) followed by a parenthesissed list of letters 
and pairs of letters. Single letters mean that symbols 
beginning with that letter are to have, the specified mode* 
unless overridden \>y a rule of higher priority, pairs, of 
letters separated by a minus sign causi? the declaration to 
apply to all lettr^riB from the first to the second inclu- 
sively. 

examples « IWPLICir REA L { I-N L INTEGER { A ~H^ 0-2) 
IMPLICIT LOGICAL (A-C,F,R.-0) 

Letters not affected by IMPLICIT stateiRents are treated 
according to the IJKLm rule. The IMPLICIT statment inay be 
placed anywhere o Its meaning inside an internal procedtire 
definition is the same as outside « 



The FUNCrXON, SUBROUTINE.* and ENTnY declarations 

The flrat two of these statements begin internal or 
external procedure definitions « The word FUNCTION or 
SUBROUTINE la followed by the nanie of the proced\xre and the 
parenthesized duiw^ symbol list. The word PUNCTION may be 
preceded by a mode word to specify the mode of the returned 
value. For an internal procedure, the word INTERN/VL precedes 
everything else» The ENTRy statement is used to specify an 
alternate entry point to a procedure ^ ENTRV is followed by 
the alternate naineo No duainy eymbol list is given , 

The END declaration 

Inside an internal procedure definition 3» the END atate- 
ment terminates the definition and re»un»s compilation of 
the main program. Otherwise, END Indicates the end of the 
aource program* 

The BLOCK mTA declaration 

The statement BLOCK lATA^ which smast be the first 
statement in the source ptoQV&m^ indicates that the program 
is a ^bXock data subprogram* « Such a program contains no 
executable statements* It consists of declarations^ parlncl"- 
pally COP^ON declarations and DATA declarations inittalizing 
variables and arrays in common storage o The block data 
subprogram is used to construct one or more common blocks to 
be linked by the loader to programs requiring them„ 



The mmhlST declaration 

Thle statement is ueed to provide format infomsition for 
10 statements that are to be processed by the nameXlst 
Interpreter » The word NftMEI*IST te followed by one or more 
declarations separated by coimnas. Each declaration is the 
name of the list enclosed In slashes followed by a list of 
variables or array names separated from each other by 
ooDsnas* This statement may be placed anywhere » Dummy symbols 
nay not appear In the llst» and the list name nsay not 
conflict with any other symbol This statement is treated In 
greater detail in the section c»i lnj3Ut/output« 

exan^le - NAMBUST /POO/A,B,C, /BAH/D,E 

The PCBMftT declaration 

This statement provides format Information for 10 state- 
ments that are to be processed by the standard format 
Interpreter* The word FORMAT Is followed by a parentheslased 
list Of format specifications and slashes <, The specif lca« 
tlons need not be separated If there 10 no ambiguity, or 
they may be separated by a comma. If too many consecutive 
xinsegarated specifications are used, the *fleld list over- 
flow'* error may occur, which can be corrected by providing 
commas. This statement miust have a statement number, and may 
be placed anywhere. It is treated In more detail in the 
section on Input/output o 

example «- 123 PCaiimT{2l4L5A3#2lA,p6.3*3c^E8.3) 



Input and Output 

An end-of-flle addreas for a READ ©tatement may b© 
specified toy placing the statement number (with no S) in 
front of the 10 list ae if it were the first element, if an 
end of file Is encountered^ control will transfer to the 
specified statement. If no end-of-flle address is specifiedf^ 
execution terminates. 

example - READ (1*123) 50, I, J, K 

will go to statement SO on end of file 

The 10 list 

The 10 list consists of one or more arrays, expressions, 
sublists, or implied DO^s separated by commas, in a WRITE 
statement the expressions may be of arbitrary complexity, in 
a HEAD statement, only variables, array elements, and araraiys 
should be used. In either case, 10 list elements must be of 
the correct mode (the mode given in the format statement). 
No mode conversion is performed. Complex data must be 
transmitted with two format specifications for real data (B, 
P, or formati)o An 10 list item should not contain a 
function call that results in the execution of another 10 
statement o 

An implied DO Is a left parenthesis, one or more 10 list 
items (including Implied DO^'s), a comnoa, index name, equals 
sign, initial value, comna, test value, comma, increment, 
and right parenthesis « If the increment and preceding comma 
are omitted, one (integer) is used. The initial value, test 
value, and increment nay be of arbitrary complexity. The 
items under control of the in^lied DO are processed repeat- 
Wly, with the index varying exactly as It would in an 
ordinary DO loop. A sublist is one or more 10 list items 
enclosed in parentheses. Implied DO^s and sublists may be 
nested to any reasonable depth. An array name is treated as 
if it were in an implied DO iterating over its linear 
subscript from beginning to end. 

exaiople - 

« PORMIVT («a table of square rootsVAoo(i5Plo.7/)) 
WRITE (0,123) (J, SftRT(PlOAT( J)), J « ^,100) 



The \inlt nvontoer 

A unit number of zero refers to the typewriter „ Other 
unit number e refer to files. 

Format etatements and arrays 

If the unit nunOder is followed by a number. It refers to 
a PCJRmT statemento If It is followed by an array name, the 
"variable format* Interpreter treats the text in the array 
as if It were a format statement. If no number or array is 
specified, the data are transmitted In linformatted "binary" 
form. The 10 list and format statement are scanned In 
parallel until the 10 list is exhausted. If the end of the 
format statement la reached, it repeats from the last left 
parenthesis « 

There are 10 types of basic forniat specifications. All 
except X and H may be preceded by a repetition count. 

Dw.d » Doubleprocl8lon« W is the field width, and d is 
the nuanber of digits after the decimal points. The 
present format interpreter ignores both. This specifi- 
cation (with its repetition count, if any) may be 
preceded by a scale factor, which is an optionally 
signed integer followed by the letter P, 

Ew.d « Heal (or one part of complex datum) « Similar to J). 
Fw.d - Real,. Similar to D. 

Ow,d - "General" (real). Similar to D. The present foroat 
interpreter does not distinguish between E, P, and 0, 

Iw - Integer (decimal) o W is the field width. The present 
format interpreter ignores it. 

Ow « Octal integer. Similar to I, 

Jm - Logical. Similar to I. 

Aw - Character. W is the number of characters. They are 
packed in an array 3 per word, 

nX - Generate n spaces in output. It is ignored in input. 
Mo data item is transmitted, 

nHstring or ^Btring« • "Hollerith", Th<& string is written 
from or read into the format statement. No data item 
is transmitted. 

A slash generates a carriage return in output. It is 
ignored in input. No data item is transmitted « 

Format specifications may be grouped within parentheses, 
A number preceding the left parenthesis will cause the 



entire group to be repeated that number of times o 

Output fomat 

Data Iteme written iinder control of I)# E» 9, 0> X« 0> and 
L fields are eeparated from each other on the same line by 
tab8« A^ Xf and H fields are not separated. L fields are 
written as t or f « 

Input f ortnat 

Leading spaces^ tabs^ and carriage returns are ignored 
when reading Dj E, P« Q» I» Oj and L data ltenis« Numbers 
mast be followed by a space « tab, or carriage return. Plus 
signs and other extraneous characters are not permitted « 
Exponential notation may not be used, h fields are scanned 
until a t or f Is founds 



Functions and eubroutines 

Functions and suloroutlnds (collectively called proce^ 
dures) are of three klnde - internal procedures^ external 
procedurea, and intrinsic functions. Internal procedures are 
those defined in the program In which they are used. An 
lntez*nal procedure definition la an arithmetic statement 
function iASF) or a section of source program beginning *flth 
the INTERNAL FUNCTIOH or INTERNAL SUBROUTINE statement and 
ending trlth the END statement. An external i»*ocedure is 
supplied toy a library, either the standard FORTRAN library 
or a library of prccediures written by the user» In the 
latter case the procedures nay be written In machine 
language or as FUNCTION or SUBROUTINE programs In FORTRAN^ 
IntarLnslc functions are simply arithmetic operators written 
m functlon-llke notation. 

Because Internal and external procedures have many fea- 
tures In common » these features will be described together. 

Function and subroutine calls 

Subroutines (internal and external) are called with the 
GPiLL statement. CALL Is followed by the name of the 
subroutine and an optional parenthesized argument list. 
Functions are called by a *f unction reference*, vrtilch Is the 
name of the function followed by a parenthesised argument 
list. In an arithmetic expression. The argument list for a 
function is required. The compiler recognizes a function 
reference by the existence of a symbol, which does not 
appear in a dimension declaration, followed by a parenthe- 
sized list. 

Function and subroutine arguments 

An argument to an internal or external procedure nay be 

(1) Any arithmetic expression. The corresponding dummy 
symbol must be of the same mode. If the expression is 
a variable or array element, the procedure may modify 
it (use it on left side of arithmetic assignment 
statement* use it in a READ statement* use it as a DO 
index, etc*), and the corresponding variable m the 
calling program will be modified. 

(S) An array name. The corresponding dummy symbol must be 
dimensioned and of the same mode. In order to make the 
subscript calculations in the procedure natch those of 
the calling progsmm, all dimensions of the dummy 
symbol except possibly the last must match those of 
the original array. For example, it is not possible to 
transmit the upper left comer of a rectangular array 
by merely using smaller dimensions in the procedure 
than in the calling program, because the first 
dimension is used for computing the correspondence 
between a given pair of subscripts and an actual 



position In the «irray. It Is poaslbXe, by using as th« 
argument an element or the array other than the first, 
to effect an offset between the array In the oalllng 
program and the duniny array In the procedure, 

(3) An external procedure name. Internal procedures and 
Intrinsic function names nay not be transmitted ae 
arguments. In tha procedure^ the dunmy symbol is used 
as the subroutine name m a CALL statement or function 
name in a function reference. In the calling program 
the external procediire must be recognized as such by 
the compiler • Use of the symbol elsewhere in the 
calling program In a CAlIi statement or function 
reference wlU suffice. If this Is not the case, the 
procedure name must appear in an EXTERNA!* declaration, 

exaoiple -^ BXl^N^L SIN 

CAUi F00(SIN,1,2) 

If the declaration were not used, the compiler would 
think SIN Is a variable, 

(4) A character string of more than 3 characters. The 
actual argument t^iat is transmitted is an array 
consisting of the characters packed 3 psi^ word, left 
justified, (A character string of 2 or fewer charac- 
ters is an integer constant and is right Justified J 

(5) A statement number. This is used for "error exits". 
It may be used only with subroutines (internal or 
external), not with functions. The argument In the 
calling program is the statement nxamber followed by 
the letter S» In the subroutine the corresponding 
dummy symbol is the times sign (x). The statement 
RETUHN n where n is an Integer and the nth dumaay 
symbol is x, will return to the statement whose number 
is the corresponding argument Instead of the statement 
following the CALL, 

Declarations on dummy symbols 

Dummy symbols in functions and subroutines must match the 
corresponding arguments. The compiler does not insert mode 
conversions for procedure arguments, and, except in the case 
of intrinsic functions, does not check for incorrect corre- 
spondence. When an argument is a function name, array name, 
or arithmetic expression, the corresponding dummy symbol 
must have the same mode, which may require a mode declara- 
tion in the procedure. If the argument Is an array name, the 
dufliny symbol must also appear in a dimension declaration. No 

v^'s^, co^MON, equivalence, namslist, internal function, 

INI^RNAL SUBROUTINE, or ENTRY declarations are permitted for 
durany symbols. The EXTERNAL declaration rasy be used where 
appllcable^k but is not necessary. 



Noym*.! return from & aiabrfxitims (tnt;«jpnal «r extern*! 5 is 
effected by the REmmM atatewant vtlth no arguatsnto Ex«cui*lo« 
reftuwes at the statement followlne the CAti:, statement;, 
Retiim to an alternate point in Wie proeraa la effects by 

RETimN n or RCTIBN (n) 

where n Is an integer and the nth dunny sy!al>oa. 1» Xo The 
aubroutlne will return to t^e atateisent whoee nimher le the 
argunent eoryeBpoaslins to the Xo 

fl 

SO ip(l ,0Eo 0) RETURN 3 

END 

If the conditlonajt in statement 20 Is satisfied^ control 
wiXI return to statement 1^3 o 

A value muet be speolfied when returning from an Internal 
or external function. This may be done either by placing the 
value in parentheses after the word aSTURN or by storing the 
value in the variable with the same najae as the function, 
and then executing a RETURN with no argument. 

example ^ INTERNAL FUNCTION POOCa,B|.C) 

RETURN Ca+B) 

POO w B+C 
RETUm^ 

END 

Inside a function {and only inside the function) the narae 
of the fwiction is treated as an ordinary variable ^ When a 
RETURN atateotent without argunsent is executed, the contents 
of that variable are used as the function value* 

Alternate entry points 

Alternate entry points to internal or external procedures 
may be specified by ui^ns of the ENTRY statement^ A 
procedure may be called hy the naine given in the Ftn^OTION or 
SUBROUTINE statement or by any nan» given in an ENTRY 
statement, in the latter case execution will begin Just 
after the ENTRY statement. 



ao TO 10 

EOTRy COS 
V"^. 57079633-1 
c ir th« varla.bXe r >i^^^ not ueed^ a«d COS storea 
c Its rsodiried argument back tnto x, the ars««^nt 
c In the calling program would be cXobber<e^ 

Ho dtuany sy»tK>l list IB given vlth %tm KJfTRY atateinent. T?he 
4ttBny ByaibcA Hat given in the TOMOriON or SUBROUTIiflS: 
atateaent applies to all entries « All &ntvy point a to a 
fwiotlon most have tho a&nim oodeo Any entry point name way 
lie uaed for atorlng the value to be rettn^ned by a fanction. 



Interrsal functions and 8ut»»mitlne« 

Internal functlona and ««broutinea ave aoapiied as p&rt 
or thm program in which th^y ar« ^sed, vf:ithout the need to 
oa IlnJced Dy the ioadero They have fche furthei' advantage 
that th«y may us& any vartabXaa in th« main program, vlfcho\«t 
cm neod for thorn to toe tocansmittad through a CQ»1MCIf block. 
They Include as a special case the •arithmetic statement 

^^*"f^i2L^. ^^^ °^'^^ ^^^^ ^^ ^ntex^isal pracedtare available in 
iioat rtmnkfi ayatems<, 

An Intex^nal pacocedure defleltidii ether than an ASF ts a 
JfJSSi??, 2L.P**^S^™ beginning with the immmt FUNCTION -©r 
lOTERNAL SUBROOTINE Stateinent and ending with the next EKB 
atatement* An explicit an^e dedar&tion cf an Internal 
fwiotlon may be made HiMedlately :following the wox^d 

examples « lOTERmL SUBROUTINE SUBH 

INTERNAI, FUNCTION FUNCa^B) 
INTERNAL LOGICAL FUNCTION LLLCj^k) 

Internal procedure deflnltione^ Including ASF definitions, 
may be placed anywhere In the source program except within 
another internal procedure definition. If an internal 
procedure definition is placed within the n»ln program, 
control will pass over It^ DO loops must be nested properly 
within each Internal procedure. It is possible to Jump out 
^Ci ^'^.i'^^^^i^J^^. procedure. J^implng into a procedure is not 
allowed. An internal procedure name may not be used as an 
argument to any procedure. 

Mode of Internal function 

The mode of the value returned by an internal function 
may be specified by a mode w<MPd in the INTERNAL FUNCTION 
statement, by an explicit mode declaration (which may be 
iS^H®^*^^ outside of the function definition h by an 
IMPLXCTT declaration, or by the IJKLMN rule. All entry 
points imist be of the same mode. 

Dummy symbols for internal procedures 

DuBKjy symbols in internal procedures are not related to 
symbols of the same name In the main program or other 
procedures. Hence a dummy ayiribol may have a type or 
dimension declaration in an Internal procedure, and the same 
ayrtbol may have a different declaration elsewhere. 

Arittoetlc statement functions 

An arithmetic statement function is a spacial type of 
internal function defined in one statement, containing the 
name of the function, the parenthesised dunwy symbol list, 
an equals sign* and an expression giving the value to be 
returned. 



An ftrlttometslo stAtement iVoncftlon is recognized by the fact 
that it appears to be an a£*lthmetlc asalgnjoent etaten^nt 
with a aubscrlpted left slde» but the ayv^iol on the left 
side does not appear in a dimension declaration. Since an 
MT definition cannot contain any declarations p the dtuongr 
BjTiii^oXa cannot he diioensionecl Ialth<High the function loay 
contain referencea to other asrrays). Ho explicit woAb 
Oeclaratlons are possible either j^ ao the modes of the d^unmy 
ayn^Xa will be deteanoKlned ^ an IKFI^ICIT declaration (if 
any)» syx> by the l^SKU^ rule. 



ExtertiaX functlona and subroutines 



Sons of the lllsrary functions (the '^standard** FORTRAN 
library functions) are known to the coffipller. This means 
that# If the name is used as an external function « the 
ooo^ler win assume the correct mode for the returned value 
unless there is an explicit mode declaration to the con- 
trmvy. No automatic mode conversion of arguoients will be 
performed. 



The standard library functl<»is are? 



arg mode value mode 



nans n 


args 


BXF 


1 


CBX? 


1 


CBXP 


1 


ALOa 


I 


DLOa 


1 


CLOO 


1 


ALCXIiO 


1 


DLOQlO 


1 


CARO 


1 


CABS 


1 


DMQD 


2 


SORT 


% 


DSQRT 


1 


CSQRT 


1 


CBitV 


1 


DCBRT 


1 


SIN 


1 


PSIN 


1 


CSIN 


1 


COS 


1 


DCOS 


1 


CCOS 


1 


TAN 


1 


DTAN 


1 


CTAN 


1 


COT 


1 


DCOT 


1 


CCOT 


1 


SINK 


1 


DSXNH 


1 


CSINH 


1 


COSH 


1 


DCOSH 


1 


CCOSH 


1 


TANK 


i 


DTANH 


1 


CTANH 


1 


ASIN 


1 


CASXN 


1 


GASZN 


1 


ACOS 


1 


mcos 


X 



description 



real real exponential 

double double double exponential 

complex complex conq^lex exponential 

real real natural logenthm 

double double double logarithm 

cosq;>lex complex con^lex logarithm 

real real common (base 10) logarithm 

double double double common (base 10} log 

coniplex real phase of conqplex number 

complex real modulus of complex number 

double double A -* [ A/b J kB 

real real square root 

double double double sqiaare root 

complex cou^lex coi^plex square root 

real real cube root 

double double double cube root 

real real sine 

double double double sine 

complex cosmlex complex sine 

real real cosine 

double double double cosine 

complex complex complex cosine 

real real tangent 

double double double tangent 

complex complex complex tangent 

real real cotangent 

double double double cotangent 

complex complex complex cotangent 

real real hyperbolic sine 

double double double hyperbolic sine 

complex complex complex lm)erbollc sine 

real real hyperbolic cosine 

double double double hyperbolic cosine 

complex coaqplex complex hyperbolic cosine 

real real hyperbolic tangent 

double double double hyperbolic tangent 

complex conrplex complex hyperbolic tangent 

real real inverse sine 

double double double Inverse sine 

complex complex complex Inverse sine 

real real inverse cosine 

double double double Inverse cosine 



C^COS i coiqplex <i(mpl^x complex inverse cosine 

ATAN i r99J. real Inverse tangant 

VK'&iim 1 double double double Inverse tangent 

CK'ffiV 1 complex complex eonplex Inverse tansent 

tiVkH& 2 real real j^iaa© of (B •♦- AlJ 

JCATANS 2 double double doubles phase of (b -h Al) 

*Jhe library hae many other functions and autorouttoea* Stor a 
complete listing of these routines j see PDP-50# RELOC^TABIiB 
SOBROOTINB LIBRARy. 

TIKI compiler anifit knoif the mode of the value of every 
external function. Bxcept In the case of ona of the standard 
library functions gSven above> a mode declaration nay be 
iiecesaary. The mode that the ooicqpller aasunes fcir the value 
of a function is determined hy explicit laode daclarationa* 
t4M known modes of standard library funotlona^ an^ impllelt 
declaratloni or the IJKU4H rule^ In tlmt order of priori ty« 

ExtexTial procedures written In FORTRAN 

An external procedure nssiy be written In FCKTHAN as a 
program whose first statement Is a FllNCTt<^ or SUBROUTDJE 
statement o The word FUNCTION or SUBROUTINE is followed by a 
parenthesized diusniy symbol list. The list is optional In the 
case of a subroutine. The laode of the value returned by an 
external function may be specified by a mode word preceding 
the word FUNCTION or by a declaration in the programo 

examples « COMPLEX FUNCTION THETA2CA,B) 

foncpjon foocx) 
sobhoutine sortCa^nI 

A FUNCTION or SUBROUTINE program ends with the END state- 
ment » as does a main program It may contain internal 
procedxares (each with Its own END statement }<. 



Xntrlyjsic functions 

An intrinsic { sometimes called '*built in**) function is an 
arithmetic operation that is used by a construct in the 
aoiirce program that looks like a function call Some 
intrinsic functions (such as ABS) are translated directly 
into machine instruct ons. Others, such as smh, are 
translated into one or more calls to library subroutines. 

The name of an intrinsic funct on may not be used in 
another context elsewhere In the program. If, for example, 
it is used as a variable, it will be a variable throughout 
the entire program, and its meaning as an intrinsic function 
will be discarded. Similarly, any dimension, OOmON» 
JSXtBSBHkhg etc, declaration, or explicit mode declaration not 
in agreement with the correct mode, will remove the meaning 
of a symbol as an intrinsic function. It is not necessary to 
declare the mode of an intrinsic function « 

Arguments to Intrinsi functions must match the require- 
ments given below for number. Automatic mode conversions of 
arguments will be performed a 



name n 


args 


arg mode 


value mode 


description 


DPLCAT 


1 


Integer 


double 


convert 


PLCAT 


1 


integer 


real 


convert 


IPIX 


i 


real 


integer 


convert 


DHU?. 


1 


real 


double 


convert 


SNQL 


1 


double 


real 


convert 


REAL 


1 


complex 


real 


real part 


Aima 


i 


complex 


real 


imaginary part 


AMOD 


2 


real 


real 


A-»[a/b]>© real remainder 


MOD 


a 


integer 


integer 


A~[A/fe]5<© remainder 


ABS 


1 


real 


real 


!A! 


lABS 


1 


integer 


integer 


|A! 


r^BS 


% 


double 


double 


|A! 


AINT 


1 


real 


real 


8ign{A)x( |A 1 ] (truncate) 


XOT 


1 


real 


Integer 


slgnCA}x|[ jA 1 ] (truncate) 


IDJNT 


X 


double 


integer 


sign(A)xt jA i j (truncate) 


SIGN 


2 


real 


real 


|A jxslsn(B) 



IS ION 


2 


Integer 


Integer 


AixsignCB) 


PSION 


a 


double 


double 


fA|xeign(B) 


DIM 


2 


real 


real 


A-min(A>B) 


TT>IM 


2 


Integer 


Integer 


A-raln(A,B) 


CMFLX 


^ 


real 


real 


A + IB 


OQ^SQ 


1 


complex 


complex 


complex conjiigate 


AWAXO 


i^ 


Integer 


real 


maximum 


kfilk'XX 


^ 


real 


real 


oaximum 


mxo 


>? 


Integer 


integer 


maximum 


mxi 


:^ 


real 


Integer 


maximum 


Dmxi 


^ 


double 


double 


maximum 


AMINO 


>? 


integer 


real 


mlnlnrum 


AMINl 


^ 


real 


real 


minimis 


MINO 


^ 


integer 


Integer 


minimum 


MINI 


^ 


real 


integer 


minimum 


DMINl 


£^ 


double 


double 


minimum 



Because of automatic mode conversion, the intrinsic 
functione to perform mode conversion (IPIX, FLOAT, SNOL^ 
DBLE, etc.) are rarely needed • In arguments to external or 
internal procedures, where automatic mode conversion does 
not take place, the use of these explicit mode conversion 
functions may be required. For example* to calculate the 
(real) square root of the Integer J, SORTCPXiOIVTCJ)) must be 
used. SQRtC^) wo\ild result in an uncorrected mode mismatch* 



APPENDIX I 

Use of the Compiler and Leader 

SlaipXe Conq;>iXatlon 

A elinple compilation ie a compilation of a naln program 
that calls no external procedures which are not stored on 
the FCHTRAN tape ($ee PDP-50, RELOCATABLE SUBROUTINE LI- 
BRARY* for a list of available subroutines and where they 
are stored). In order to compile a FORTRAN main program It 
is necessary to have the microtape containing the FORTRAN 
compiler mounted on tape drive n» The compilation is begun 
by typing the comand 

nP 

to ID or to ET (For complete descriptions of these programs, 
see PDP-23* INVISIBLE DEBUCKIER CXD), and PDP-S£S* EXPENSIVE 
TYPEWRITER Jet))o The compiler will type the title (first 
nonempty line) of the prog3?am on the typewriter as it begins 
the compilation. Any errors detected by the compiler will 
cause informative messages to be printed » 

If the oompllation falls, control will be returned to ID. 
If the conogpilation is successful^i the loader will be started 
automatically. When the loader types "options « "^ type a 
carriage return. If loading ie successful, the siae of the 
loaded program will be printed. Otherwise, an error message 
will be printed. The loaded program will be placed onto drum 
field 1 and into core. If the loaded program is more than 
4096 words long, only the first 4096 words will be placed 
onto drum field ±, Control is passed to ID. 

Execution of the program nay be begun by typing 



or 

ioea 

to ID, If the program crashes, control will be returned to 
ID which will type a message of the form? 

xxxjrxtt hit (These are ID*^s standard messages for in- 
or indicating that a halt instruction or an 

yyyyy« zzzzz illegal memory reference has occurred) 

To stop a running program* hit the call button. This will 
also return control to ID. From ID, control may be returned 
to ET by typing 

£ 

or a fresh copy of the program (or the first 4096 words of 



th© j^rogram if the program 1« longer than 4096 words) laay be 
obtained and execution started by typing? 

(zero underl)ar) 

IfflU 

lOlSQ 

Coa^lling a Subroutine or a Function Subprogram 

The proced\ire for compiling a subroutine or a function is 
the aame as that for conqpiling a main program, except that 
at the termirmtion of compilation tho loader is not started 
but rather control is passed to ID. If the compilation was 
successful^ the object file can be placed in a library* IPhls 
is dcane by saving the object file on any microtape contain- 
ing a file system (any of the public tapes numbered % 
through 7 will do)<. The object program can be saved by 
typing the following sequence of coiaraands^ 

nF (n is th® number of the drive on which the file 

tape is mounted) 
8 fortllb progname Cprosr«ame is any desired file name 

(nine or fewer characters)) 
b 

(For an explanation of what these commands do^ see PDP-42, 
MICROTAPE PILE SYSTEM)^ At the end of this sequence ID is in 
control and the compiled subprogram is stored on the 
microtape mounted on drive n« The next section explains how 
to compile a main program which calls subprograms saved on a 
library tape^ 

General Compilation 

To compile a main program which uses subprograms that are 
saved on a library tape 5 proceed as in the case of a siirq;)le 
compilation* At the end of a successful compilation^ the 
loader program will be started* When the loader types 
"options ~ "^ type the number of the tape drive on which the 
libary tape is i?K>unted folicswed by a carriage return* The 
loader will then search all files on that library tape whose 
first name is ^fortllb" for procedux'es that have been 
referenced „ One external procedure n»y reference other 
external procedures o All will then jproceed as in a simple 
compilation o 

lioader 

The loader combines relocatable files into binary pro- 
grams. The operation of the loader may be tmire closely 
observed by the use of various options which may be selected 
when the loader types ^options • *, At this pointy various 
characters may be typed* Each character specifies an option* 
If a character \^ich does not correspond to any option is 
typed, the loader will type "options • ** again* 



CHARACTER OPTION 

any digit The digit Is taken to t>e the \init niwiber of 

a tape drive containing a library tape« 

m Print a load map giving the name-, location, 
and length of each procedure that is 
loaded. 

r Print every reference to an external proce* 
dxareo This option la only effective if 
the load map option is also selected. 



Features Not lmplemente<l 
Non-6tandard sutjroutlna returns {i*e, RETURN 3), 
BLOCK CATA subroutines , 
NAMELXST X0» 
VftrlttbXe fomftt 10, 
Entlra arrays In 10 listo 
Conqpldx item in XO Iiat« 

10 to devices other than the typewriter through WRrTE and 
HEAD statements » 

Scale factors (P field } in FORmT Btatesiente. 
H format input. 

Use of an internai function name as a "variable* inside of 
the Internal function definitions in any context other 
than the following 2 

NA^E » expression 
RETUHIN 

The following intrinsic functions are not implemented? 

DIM 

torn 

Amxi 
mxi 

AWXKl 
MINI 
DMAX^ 
DMIHl 



Known Biss® 

doablep^^eclslon t real does not work. The neciss^jary aonmr^- 
«ton falls to occur. 

The usQ of Interml r^jiBctions In expr«ssloBs is dangeroua. 
T^mgoT^vX^^ used toy the staterfiaot calling the iz^teraal 
fttooticn my J)^ altered aut'li^g the avaluatlon of fches 
Internal furscfclcn. 

Exceedingly long cnaractsr ^string constants (l.a. thousaiids 
of characters) vill net oompll^, 

n^mhH coded extpr^rri&l p2*oc?edure8 {functions csr subroutines) 
wilx not tfork correctly if th&y are loaded tnto core 
modules pfcher than o &.mJ they contain arrays that are 
nelr.her in common nor parameters to the prooediire. th© 
loaa nap option of the loader can be \m^d to detect thla 
coi^ltion. mtl^ hus may ha avoided toy putting all arrays 
occuring In FOI^BAN coded subprograms in QOWIQ}^ blocks. 



